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^ ■ ABSTRACT 

I , We report the results of a spectropolarimetric survey of a complete far infrared se- 

^ ■ lected sample of Seyfert 2 galaxies. We have found polarized broad Ha emission in 

' one new source, NGC5995. In the sample as a whole, there is a clear tendency for 

, galaxies in which we have detected broad Ha in polarized light to have warm mid-far 

' infrared colours (Fgo^ni/F25/im ^4), in agreement with our previous results. However, 

a comparison of the optical, radio and hard x-ray properties of these systems leads us 

, to conclude that this is a secondary consequence of the true mechanism governing our 

' ■ ability to see scattered light from the broad line region. We find a strong trend for 

I galaxies showing such emission to lie above a critical value of the relative luminosity 

f~| . of the active core to the host galaxy (as measured from the [OHI] 5007A equivalent 

C^h' width) which varies as a function of the obscuring column density as measured from 

Q , hard x-ray observations. The warmth of the infrared colours is then largely due to a 

■ combination of the luminosity of the active core, the obscuring column and the relative 

' importance of the host galaxy in powering the far infrared emission, and not solely 

^ , orientation as we inferred in our previous paper. Our data may also provide an expla- 

IL" ' nation as to why the most highly polarized galaxies, which appear to have tori that 

. 5^ I are largely edge-on, are also the most luminous and have the most easily detectable 

^ . scattered broad Ha. 

^ . Key words: galaxies: Seyfert - galaxies: active - polarization - scattering - infrared: 

■ - - ' galaxies - X-rays: galaxies 



1 INTRODUCTION 

In the standard model for active galaxies (hereafter AGN), 
hard continuum radiation from the accretion disk around a 
massive black hole ionises the surrounding gas, which can be 
divided into two components: the broad line region (BLR) 
and the narrow line region (NLR) , at distances from the core 
of < Ipc and < Ikpc, respectively (cf Osterbrock 1993). For 
the case of Seyfert galaxies, the classification into Seyfert 1 
or 2 depends solely on the presence or absence of consider- 
ably broadened (several thousand kms~^) permitted lines: 
ie, whether or not we can see the BLR in terms of the stan- 



* Charlene Heisler died after a long struggle with ill-health on 
October 28 1999. She will be greatly missed by all who knew her. 



dard model. The simplest form of 'unification' for Seyfert 
galaxies invokes a dusty torus to obscure the light from the 
BLR in the Seyfert 2s, with the sole free parameter being 
inclination (see, for example, Antonucci 1993 for a review of 
unification schemes for Seyfert galaxies). 

The firmest evidence in favour of this simple unified 
model comes from optical spectropolarimetry and x-ray 
spectroscopy. Spectropolarimetry reveals scattered broad 
permitted lines (eg. Antonucci and Miller 1985, Miller and 
Goodrich 1990). This can only come from a 'hidden' BLR 
(hereafter HBLR), that is obscured from our direct line of 
sight. Similarly, x-ray data for Seyfert 2s show evidence for 
either an absorbed Seyfert 1 like continuum in the hard x-ray 
band, or when the extinction is sufficiently large, a scattered 
spectrum characterised by an Fe Kq line with large equiv- 
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alent width (sec, for example, the recent summary of the 
results from the ASCA satellite by Turner et al. 1997a,b). 
The evidence that at least some Scyfcrt 2s are in fact 'mis- 
aligned' Seyfert Is is therefore compelling. However, previ- 
ous samples in both the optical and x-ray have been selected 
in a fairly ad-hoc fashion. For the x-ray the limitation was 
simply sensitivity, so only the brightest x-ray sources, which 
also tend to be the least obscured, were studied. For the op- 
tical spectropolarimetry, sources were often initially chosen 
on the basis of their known high continuum polarization (eg. 
Miller and Goodrich 1990) or from ad hoc samples drawn 
from the literature (eg. Tran, Miller and Kay 1992, Kay 
1994, Young et al. 1996a). In addition, with the exception 
of Kay (1994), who did not cover Ha in the spectral range 
of her data somewhat limiting its usefulness, and Young et 
al. (1996a), non-detections of HBLRs were rarely published. 

Therefore it is not possible to use the published results 
on optical spectropolarimetry to test the unified model in a 
meaningful statistical sense, or, indeed, to determine if the 
simplest form of unification is actually the best match. In 
order to interpret whether the detection of an HBLR relates 
to orientation, or perhaps to other properties of the host 
galaxies, we must first understand how the irutial galaxy 
sample is chosen. We therefore undertook a new survey in 
which we obtained optical spectropolarimetry of a sample 
with well understood far infrared properties. As discussed 
in Section 2, the far infrared luminosity is likely to be a rel- 
atively extinction (and hence orientation) independent mea- 
sure of the combined AGN and host gala^xy luminosity. In 
addition, there was a notable trend shown in early detec- 
tions of HBLRs for those galaxies to have 'warm' mid-far 
infrared colours (eg Inglis et al. 1993). We therefore wished 
to determine if this was still true for a larger sample, and 
that the AGN with cooler colours did not show evidence 
for HBLRs. The initial results of that survey were reported 
in Heisler, Lumsden and Bailey (1997: hereafter Paper 1). 
There we reported an apparent trend for HBLRs to only 
appear in those AGN where F60pm/F25;jm < 4. Since the 
mid-fax infrared emission in Seyfert galaxies is dominated 
by thermal emission from dust, this trend is consistent with 
the unified model, since low values of the Ffi()^m/F25pm ratio 
imply warmer dust temperatures. Therefore we concluded 
that galaxies with warm mid-far infrared colours had tori 
that were closer to face-on, and hence more likely to show 
scattered radiation from the BLR. 

However, our initial survey did suffer from some prob- 
lems relating to the actual sample observed (see Section 2 for 
details), and was simply too small to test for other depen- 
dencies in the data that might mimic the effect that orien- 
tation would have. In particular, it is possible that star for- 
mation may play an increasingly important role in the AGN 
with cooler mid-far infrared colours, since the dust emis- 
sion surrounding HII regions typically peaks nearer 100/xm. 
Alexander (2001) suggested that the hard x-ray emission 
from the sample given in Paper 1 is consistent with such 
a picture, since the AGN with both warm and cool colours 
show a similar spread in implied neutral hydrogen column 
density. He inferred from this that the mid-far infrared 
colours related more to the presence or absence of signifi- 
cant star formation rather than orientation. In addition he 
noted, as had Tran et al. (1999), that there was an apparent 
trend for the HBLRs to have higher [OIII] 5007A/H/3 ratios, 



and indeed for this ratio to be correlated with mid-far in- 
frared colour. Again, this may be consistent with enhanced 
star formation in the galaxies with cooler colours, since the 
starburst will enhance H/3 but have relatively little effect on 
[OIII] emission in most cases. 

This paper therefore presents the results of an extended 
survey of fax infrared selected AGN. Since Paper 1 we have 
also refined our selection to exclude objects whose actual 
classification is in doubt (see Section 2.5). We will present 
the full spcctropolarimetric data, including non-detections, 
for all objects for which we acquired data in a later paper. 
We will also present detailed modelling of all aspects of the 
observed polarimctry rather than just a simple considera- 
tion of the presence or absence of scattered broadened Ha 
emission in that paper. 

2 THE SAMPLE AND OBSERVING DETAILS 
2.1 Selection Method 

Any test of AGN unification must use a sample selected on 
a known isotropic property since the unified model implies 
an orientation dependence for the type of Seyfert seen. It 
is important to avoid the error, often seen in the literature, 
that all Seyfert 2's are identical. If the unified model is cor- 
rect, Seyfert 2s must have a wide range of orientation of the 
AGN axis with respect to our line of sight, from narrowly 
mis-aligned so that the Seyfert 1 core is only just out of 
view, to those where the broad and narrow line regions are 
almost entirely in the plane of the sky. It is highly unlikely 
that the observed properties of Seyfert 2's seen at widely 
different orientation angles are the same. Therefore it is not 
sufficient to pick a sample of Seyfert 2 galaxies at random 
from the literature, as any such sample will almost certainly 
be biased with respect to orientation. 

In practice, a combination of extinction, opacity effects, 
and varying beam-size at different wavelengths, makes it im- 
practical to select on any property related solely to the ac- 
tive nucleus. It is possible though to select based on intrinsic 
properties of the galaxy as a whole, in a wavelength regime 
where the active rmcleus itself is likely to be important. Ex- 
amples of this type include far infrared (hereafter FIR), cm 
radio, and global optical emission. Each of these has its own 
particular advantages and disadvantages, but we initially 
chose a FIR selection since there is a complete FIR selected 
galajcy database readily available. 

The nature of the FIR emission in Seyfert galaxies was 
the subject of considerable speculation until relatively re- 
cently. The turnover in the spectral energy distribution be- 
yond 100/um was seen as evidence for a thermal origin (Chini, 
Kreysa and Biermanu 1989). More recent ISO photometry 
has helped to confirm that this picture is correct for most 
if not all nearby radio quiet AGN (eg Perez Garcia and Ro- 
driguez Espinosa 2000). There have been many attempts to 
model the properties of the dust emission around the ac- 
tive nucleus itself in the context of the unified model (eg 
Pier and Krolik 1992, Laor and Draine 1993, Granato and 
Danese 1994, Efstathiou and Rowan- Robinson 1995). These 
models suggest that the 60-lOO^m emission from any AGN 
is largely unaffected by orientation. 

Some fraction of the emission at these longer wave- 
lengths could be due to the host galaxy rather than the core 
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however. Star formation in particular gives rise to significant 
FIR emission. It is possible however to distinguish the main 
source of the FIR emission using simple colour-colour dia- 
grams of the observed IRAS fluxes (eg. Dopita et al 1998). 
Figure 1 shows IRAS colour-colour plots for the sample of 
galaxies considered by Kewley et al. (2000). That sample 
shares similar FIR selection criteria to those used in Paper 
1 and this paper, but included all types of galaxy. The data 
shown here differ from that in Kewley et al. since we have re- 
vised the IRAS fluxes to use the data from the IRAS Bright 
Galaxy Survey (BGS: Soifer et al. 1989 and Sanders et al. 
1995). The solid line in the figures represents the track for 
an increasingly reddened Scyfcrt galaxy, whereas the dashed 
lino shows the locus of known starburst galaxies (see Dopita 
et al. 1998). Galaxies of mixed excitation should lie between 
these lines (the dot-dashed line shows the extreme mixing 
line between starburst and AGN activity). Seyfert galaxies 
of both types are shown in the left hand panel, and star- 
bursts and LINERs in the right, in order to show the clear 
differences in their observed properties. 

2.2 The Infrared Catalogue 

We used the catalogue of IRAS galaxies and the redshift 
survey of Strauss et al. (1990, 1992) to identify our AGN 
sample. We did not rely on the Rush et al. (1993) cata- 
logue used in Paper 1 for several reasons. First, it suffers 
badly from classification errors. This meant we had to in- 
dependently classify all galaxies that met our selection cri- 
teria regardless of which source we used for the IRAS data. 
Secondly, the Strauss et al. catalogue is larger, and com- 
plete to fainter flux limits, making it better for identifying 
a larger parent sample of AGN. Lastly, doubts have been 
raised about the accuracy of the IRAS fluxes in the Rush et 
al. catalogue (Alexander & Aussel 2000). 

Strauss et al. (1992) carried out a spectroscopic survey 
of sources selected from the IRAS Point Source Catalogue 
(PSC v2.0) by Strauss et al. (1990), to determine redshifts 
of the galaxies, and remove non-extragalactic sources. The 
actual IRAS fluxes they quote are not drawn solely from the 
PSC however. Where the PSC flags a source as extended or 
variable, Strauss et al. (1990) derived ADDSCAN fluxes. 
These should in principle be better than the PSC fluxes 
since (i) they correctly sum all the flux from an extended 
object and (ii) they coadd more of the raw IRAS data and 
are therefore of higher signal-to-noise. Strauss et al. (1990) 
discuss this issue at greater length. The final catalogue of 
Strauss et al. (1992) contains most galaxies in the original 
PSC down to a 60/im flux limit of 1.9Jy. 

The IRAS fluxes actually quoted in this paper however 
are derived from a variety of sources, since we sought the 
best possible data once the initial selection was carried out. 
We compared the fluxes in Strauss et al., the IRAS Faint 
Source Catalogue Version 2 (FSC: Moshir et al. 1991), and 
the BGS (Soifer et al. 1989, Sanders et al. 1995). These 
all agree within the quoted errors for point sources. The 
BGS data are also ADDSCAN fluxes of the original raw 
data, but with a higher degree of rejection of 'bad' scans. 
They should therefore be more accurate than the Strauss et 
al. data. Therefore, for galaxies with Feo^jm > 5Jy we used 
the fluxes from the BGS. We used the data of Strauss et 
al. (1990) for galaxies with Feo^m < 5Jy where those data 



were derived from ADDSCAN fluxes. Finally, for all other 
galaxies we used the data quoted in the FSC (which again 
is produced from coadding more of the raw IRAS data than 
the PSC). The BGS data quote errors on their fluxes, as 
does the FSC. Wo assume the errors on the Strauss et al. 
ADDSCAN fluxes must be bound by the FSC errors. 

Our initial sample is based on those galaxies with good 
detections in the PSC (ie quality flag set to 3) at 60/im. They 
also have moderate or good detections at 25 and lOO/nm 
(quality flag set to either 2 or 3). Wc did not set any re- 
striction on the 12fim data. However, because we do not use 
the actual PSC data in our analysis, all objects in our fi- 
nal sample actually have firm detections at all four IRAS 
wavebands. 



2.3 Detailed Selection Parameters 

The detailed selection criteria of our sample are as fol- 
lows (some of these arc actually imposed by using the 
Strauss et al. catalogue): Lfir > 10"'^'' Lq( where Ffir has 
its usual definition - Ffir = 1-26 (2.58F60Mm + FiooMm) x 
10~^'*Wm~^); F60nm/F25,:jm < 8.5 which matches most 
known Seyferts (de Grijp et al. 1992) without requiring us to 
classify the large population of normal galaxies with cooler 
IRAS colours; Feo^m > 3Jy; Galactic latitude |fe| > 25°; 
declination less than -|-70°; and classification as a Seyfert 
2 (from our own spectroscopy if possible: see Section 2.5). 
We note that we class Seyfert 1.8 and 1.9 galaxies as 
Seyfert 2's for the purpose of this paper. We use a value 
of Ho — 75kms^^Mpc^^ throughout. 

These criteria are looser than those adopted in Paper 
1. Despite this, some of the galaxies discussed there do not 
fall into the current sample. These are NGC34 (more accu- 
rately classified as a LINER), NGC7496 and NGC7590 (both 
Lfir < W^^ L©)- None of these objects showed evidence for 
an HBLR. 

Table 1 lists the final sample of 28 galaxies that satisfy 
all our criteria. We do not have spectropolarimetry for 4 of 
these objects. All have cool colours, with Ffio^jm/F^r.jnu > 6. 
A detailed discussion of these sources is deferred until the 
Appendix, though we will note, where appropriate, the effect 
they have on our analysis. 



2.4 Observations 

Full details of the observing procedures and data reduction 
will be given in a future paper. With the exception of those 
objects for which wc have used data available in the litera- 
ture, all spectropolarimetry was obtained at either the AAT 
or WHT, on the nights of 16, 17 and 18 August 1995 (AAT), 
29 July 1996 (WHT), 26, 27 and 28 August 1997 (AAT) and 
7, 8, 9 and 10 March 1998 (WHT). Only the data from 1998 
were obtained with the slit at the parallactic angle (mainly 
for reasons of practical convenience when carrying out spec- 
tropolarimetry). We centred the slit on the brightest part 
of the galaxy at approximately the V band. The conditions 
were photometric for the all except the night of 18 August 
1995. Total exposure times were approximately two hours 
for all sources observed at the AAT, and one hour for the 
sources observed with the more efficient ISIS spectrograph 
at the WHT. Ideally, the observations would seek to achieve 
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a level consistent with that expected of any broad line in 
the scattered light (cf Section 3.1). In practice, instrumen- 
tal and other effects limited us to a fixed small error in the 
polarization of ~ 0.3%. The fainter sources, and some that 
were observed in adverse weather conditions, have larger ob- 
served errors. Fuller details of the errors will be given in the 
paper presenting the full data. 

For the AAT data, we used the RGO Spectrograph to- 
gether with a Tektronix 1024^ CCD. For the WHT data, we 
used the red arm of the double beam spectrograph ISIS with 
a Tektronix 1024^ CCD. A calcite prism was used in both 
instruments to split the incoming beam into e and orays to- 
gether with an aperture mask made up of discrete slit-lets. 
A half-waveplate modulated the incoming phase. Four steps 
of the waveplate at 0°, 45°, 22.5° and 67.5° were required to 
derive the full set of Stokes parameters for linearly polarized 
light. The slit width was typically 2 arcseconds at the AAT, 
and 1-1.5 arcseconds at the WHT. We extracted data from 
a region of between 3 and 4 arcseconds along the slit in most 
cases. Sky subtraction was achieved by nodding the object 
into an adjacent slit-let. The data were reduced in a stan- 
dard fashion. The effective spectral resolution is between 7 
and 9A for all data depending on the exact slit width used. 
The sky and object spectra for the four separate waveplate 
positions were extracted then combined to give the final Q, 
U and I Stokes parameters. 

2.5 Spectral Cleissification 

We measured the relative line fluxes from our spectra by fit- 
ting Gaussian profiles. Where the actual observed profile is 
highly non-Gaussian, we used a multiple component Gaus- 
sian fit to measure the total flux, since we are not interested 
in the detailed line profile information itself. We also allowed 
for the possibility of stellar H/3 absorption in our spectra. It 
is possible to overestimate the internal extinction, and mis- 
classify starburst galaxies as AGN without this correction. 
We felt it was better to actually fit the absorption compo- 
nent directly rather than making an average correction as 
has often been done in the past. This will make little differ- 
ence for galaxies with an old underlying stellar population, 
where the H/3 absorption is weak in any event, but has a 
significant effect on those galaxies which have a stellar con- 
tinuum similar to that of an A star (essentially those galaxies 
with a post-starburst population). 

The extinction was derived using the Whitford redden- 
ing curve as parameterized by Miller & Mathews (1972), 
assuming the intrinsic Ha/H/? ratio is 3.1. The extinction 
corrected line ratios for our own spectra, and for those galax- 
ies taken from the literature, arc given in Table 2. We note 
that the derived E{B — V) values may not be completely 
accurate since we did not observe at the parallactic angle 
for most of our observations. This is especially a problem 
for strongly nucleated sources observed well away from the 
parallactic angle, since a point source at the wavelength of 
H/3 will appear to shift by up to 1 arcsecond from its position 
at Ha due to differential atmospheric dispersion. Similarly, 
such effects can lead to errors in line fluxes that have been 
corrected for extinction. By contrast, line ratios, as long as 
they arise from the same component, are largely unaffected, 
since the extinction correction essentially forces the data to 
have the 'correct' spectral form. The same arguments can 



also be made regarding data taken under non-photometric 
conditions. These caveats also apply to much of the data 
presented in the literature, and should be borne in mind as 
a possible source of error in the extinction corrected [QUI] 
5007A line luminosities quoted in Table 2 and used through- 
out this paper. 

We classified each galaxy based on the diagnostic dia^ 
grams of Veilleux and Osterbrock (1987). We also considered 
whether the galaxy showed evidence of Hell 4686A emis- 
sion. The high signal to noise of our spcctropolarimctry data 
means we are able to detect this line where it is present. This 
line is, by definition, not seen in LINER galaxies. It is seen 
in some starburst galaxies, where it arises in the winds of 
Wolf-Rayet stars (see, eg, Vacca and Conti 1992). However, 
the feature in Wolf-Rayet galaxies is much broader than the 
other observed lines, whereas in all the galaxies we observed 
the width of the 4686A line is consistent with the width of 
[QUI] 5OO7A, indicating it also arises in the NLR. 

Four of the galaxies included in Table 2 do not have 
standard Veilleux & Osterbrock Seyfert classification. We 
have kept these galaxies in our sample on the basis of other 
reported characteristics. 0031-2142 is classed as a Seyfert 
1.8 by Moran et al. (1996) on the basis of weak broad Ha 
emission. Mrk 334 is classed as a Seyfert on the basis of a 
weak [FeX] fine, which cannot arise in a LINER. 0019-7926 
shows weak Hell emission. It also seems possible that we 
have not observed the actual nucleus. Both de Grijp et al. 
(1992) and Vadcr et al. (1993) present spectra of this galaxy 
which have a standard Seyfert classification on the basis of 
all the Veilleux & Osterbrock (1987) indices, as well as signif- 
icantly stronger [OIII] 5007A emission than is present in our 
spectrum. However, Kewley et al. (2001) present data that 
arc in agreement with ours. Given the uncertainty we have 
decided to include this object. NGC 7582 shows weak Hell 
emission, and has been observed to show transient broad 
Ho (Arcxtaga ct al. 1999). We tested for any effect that the 
inclusion of these galaxies had on the statistical tests that 
we describe in Section 3. There was no significant difference 
in any case. The only difference that their exclusion makes 
is to the overall detection rate reported in Section 3.1, and 
we discuss that factor there. We also discuss, where appro- 
priate, any other effect that the inclusion of these galaxies 
may have. 



3 RESULTS 

3.1 The Polarimetry Data 

Only one of the galaxies newly reported in this paper, 
NGC5995, contains an obvious HBLR. Broad Ha and Hp 
are clearly evident in the polarized fiux from this object 
(see Figure 2). Here we have plotted the Stokes flux from 
our observations (derived from the rotated Stokes parame- 
ter - these quantities are a less biased estimate of the actual 
polarization and polarized flux present since they are not 
positive biased - see Miller, Robinson and Goodrich 1988). 
The measured full width at half maximum of the broad Ha 
component is approximately 2700kms~^. Although the data 
has lower signal-to-noise near H/3 we can also measure the 
Ha/H/3 ratio in our polarized flux spectrum. We obtain a 
ratio of 18±6, which corresponds to E{B - V) = 1.7 ± 0.4 
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assuming an intrinsic ratio of 3.1. We can compare this 
with the extinction value given in Table 2, which shows 
E{B — V) = 1.78. Given the similarity of these values it 
would appear that most of the extinction lies between us 
and the scattering particles or narrow line region, rather 
than between the AGN core and the scatterers. Finally, one 
other notable feature of the spectrum of NGC5995 is that it 
shows clear evidence for strong stellar absorption lines. This 
is rather unusual for HBLRs as noted by Kay and Moran 
(1998). 

Two other galaxies show weak features that may be the 
signature of an HBLR. NGC5135 has a weak broad feature 
in the polarized flux, but no signature in the polarization 
itself. In addition, the residual broad feature left after sub- 
tracting out the narrow hue component is considerably red- 
shifted with respect to the systemic velocity (by more than 
lOOOkms"^). This could be due to the weakness of the fea- 
ture however, since it is possible that wc have ovor-corroctod 
for the narrow Ha and [Nil] components in the polarized 
flux spectra, and hence removed the blue 'wing' of the broad 
feature. NGC5929 shows a weak broad feature in the polar- 
ization at Ha but no evidence for a broad component in the 
actual polarized flux. Again, it is possible that a weak broad 
line could be masked by the much stronger polarized narrow 
lines. We do not therefore consider these features as detec- 
tions of HBLRs in either galaxy, but clearly both sources 
are worthy of further study to confirm whether this is true. 

We also considered whether we would expect to see an 
HBLR given the signal-to-noise achieved in these observa- 
tions. We derived predicted broad Ha fluxes from the ob- 
served [OIII] 5007A fluxes, with due correction for the dif- 
ference in extinction between the two lines, but assuming 
that both the [0111] emission from the NLR and the scat- 
tered light from the BLR suffer the same extinction along 
our line of sight. We used two methods for deriving the pre- 
dicted flux. We fitted the observed relation between [OIII] 
5007A and broad Ha for the Seyfert Is observed by Stirpe 
(1990), and then assumed that an average of 2% of the BLR 
light is scattered into our direction in Seyfert 2s (see Lums- 
den & Alexander in preparation) for the first prediction. For 
the second we fitted the observed relation between the scat- 
tered broad Ha flux and the direct [OIII] 5007A flux in the 
HBLRs studied by Tran (1995a), Young et al. (1996a) and 
ourselves. The two results agree to within a factor of 5. We 
then derived an upper limit to the observed scattered broad 
Ha in our data from the Stokes flux, assuming the full width 
at half maximum of the line was 3600kms~^ (the average for 
the Seyfert Is in Stirpe 1990). Where our limit is 3 times 
less than the predicted value in both cases we assume that 
the non-detection is secure. Where this condition is satisfied 
by only one of the predicted fiuxes, we consider the non- 
detection is likely but flag the possibility that this is wrong 
in Table 2. In two cases, NGC7479 and 1408-1-1347, our hmit 
is above the predicted values, and we cannot therefore rule 
out the presence of an HBLR. These results are broadly in 
line with the analysis of Alexander (2001) of the sample in 
Paper 1. Further observations of all the galaxies for which 
we cannot absolutely rule out the presence of an HBLR are 
clearly desirable. For the rest of this paper, however, wo will 
assume that all of the galaxies in which we did not detect 
an HBLR are in fact non-HBLRs. 

In total 8 out of the sample of 24 observed galaxies show 



evidence for a HBLR. Assuming the 4 unobserved galax- 
ies with cooler IRAS colours arc not HBLRs, this gives an 
overall detection rate of 29%. If we exclude the 4 galaxies 
with predominantly intermediate or LINER types discussed 
in section 2.5, the detection rate would be 8 out of 21 ob- 
served, and a total sample of 24. This gives a rate of 33%. 
Both of these are close to the value reported by Moran et al. 
(2000) based on an initial volume limited optical selection of 
35%. The similarity probably indicates our initial selection 
is largely unbiased and provides a fair representation of the 
full range of nearby Seyfert 2s. 

3.2 Optical and Infrared Diagnostics 

Figure 3 shows where the HBLRs in our sample fall as a 
function of IRAS colour. Our expanded survey indicates that 
the trends we reported in Paper 1 still hold. There is a clear 
correlation between the mid-far infrared flux ratio and the 
ability to detect an HBLR. Galaxies with warmer colours 
(F60|tim/F25^m < 4) almost uniformly show evidence for an 
HBLR: the two exceptions are the optically faint galaxies 
0019-7926 and 0425-0440. As noted in Section 2.5 we may 
not have observed the nucleus of 0019-7926. If we have in- 
deed misclassified this particular galaxy the net effect on our 
results is small as discussed in more detail in the Appendix. 
0425-0440 is one of the galaxies in which we cannot rule 
out the possibility of a misclassification due to insufficient 
signal-to-noise. 

Figure 3 can also be compared directly with Figure 1. It 
is clear that there is a greater spread in the Seyfert sample 
under consideration here than in the Kowloy et al. (2000) 
sample. This is largely due to the inclusion of the galaxies 
which do not have standard Seyfert classifications (see Sec- 
tion 2.5), the low luminosity Seyferts in NGC5194 (M51) and 
NGC7172, and the inclusion of the galaxies for which we do 
not have spectropolarimetry (for which wc arc generally re- 
liant on only one literature source for classification). In most 
of these cases it is likely from the IRAS data alone that the 
active core is not the dominant source of luminosity, since 
they follow the same track as the starburst and LINER pop- 
ulation in Figure 1. By contrast the HBLRs uniformly follow 
the reddened Seyfert 1 track. The other striking aspect of 
the HBLRs from Table 2 is that they are always classified as 
Seyferts on the three standard Veilleux & Osterbrock clas- 
sification diagrams. The HBLRs and the starburst/LINER 
dominated galaxies represent the two extremes in our data. 
The rest of the non-HBLR population lies between the two, 
with some on the Seyfert 1 reddening track, and some on 
the starburst/LINER track. 

We considered whether luminosity is a factor in the de- 
tection of the HBLRs. A Mann- Whitney U test shows that 
the distributions of Lfir for the HBLRs and non-HBLRs 
are consistent at the 95% level. Therefore Lfir is not a fac- 
tor in determining the presence of a HBLR, in agreement 
with our findings in Paper 1. A different result is obtained if 
we look at the [OIII] 5007A line luminosity, Lasoo?. It is gen- 
erally held that the luminosity in this line is a good measure 
of the luminosity of the active core (Mulchaey et al. 1994; 
Alonso-Hcrrcro, Ward and Kotilaincn 1997). Our data sug- 
gest marginal support for there being a difference between 
the HBLRs and non-HBLRs at the 10% significance level, 
in the sense that the HBLRs are more luminous on average. 
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If we look at how the IRAS colours behave as a function 
of La5oo7, wc find that all galaxies with Lasoot > 10® L© lie 
near the Seyfert 1 reddening line in Figure 3, regardless of 
whether they are HBLRs. This last point is confirmed by the 
weak correlation between Lasoot and IRAS colour at a 5% 
significance level shown in Figure 4(a): the more luminous 
sources have warmer colours. There is no such correlation 
between Lf/h and colour as shown by Figure 4(b). 

Another way to demonstrate the dependence of the 
IRAS fluxes on the AGN luminosity is to consider what frac- 
tion of the mid infrared flux actually arises near the active 
core. We have used the small aperture 10/xm data presented 
in Giuricin, Mardirossian and Mozzctti (1995) and Maiolino 
et al. (1995) to determine how 'compact' the IRAS 12/im 
emission is. This is done in the standard fashion by deflning 
a compactness parameter, C, which is the ratio of the small 
aperture lO/um flux and the IRAS 12/xm flux scaled by an 
appropriate colour correction (Devereux 1987). The results 
are given in Table 3. Values of C near 1 indicate that the 
IRAS emission arises from within the aperture of the lO/itm 
beam. Although only half of our sources have small aperture 
10/xm data, the results are instructive. All of the galaxies 
with La5oo7 > 10* '^ L0 have C> 0.65. Only NGC5135 has a 
value of La5007 < 10*'® Lq and a value of C larger than 0.65. 
These results are consistent with our findings above, since 
they indicate that the galaxies with lower AGN luminosities 
are more likely to show significant host galaxy emission in 
the mid infrared, which must also be true in the FIR as well. 

We also used the equivalent width of the [OIII] line, 
Wa5007, as a measure of the ratio of the luminosities in the 
active core and the host galaxy. Figure 5 shows that Wa5007 
and La5007 are correlated as might be expected, but only 
at the 10% significance level, and there is a large scatter. 
Galaxies with high values of Wa5007 tend to be more lumi- 
nous, and to be more clearly classified as Seyferts (Table 2). 
In this sense, Wa5007 is a good indicator of whether or not 
the AGN dominates the optical spectrum and the overall 
spectral energy distribution. 

In Figure 6(a) we show the relation between Wa5007 
and F60Mm/F25^m. Those galaxies which have large equiva- 
lent widths also show HBLRs. Formally, the galaxies with 
and without HBLRs have distributions of Wa5007 that are 
different at the 99% confidence level. There axe counter ex- 
amples, since NGC5995 and 0618-2524 both have values 
of Wa5007 which are more typical of the bulk of our sam- 
ple. Those galaxies with the largest Wa5()()7 also have the 
warmest colours. Again this is consistent with the picture 
outlined above with regard to AGN luminosity and mid-far 
infrared colour. 

As noted by Tran et al. (1999) and Alexander (2001) 
there is a clear trend for Seyfert 2s with HBLRs to have both 
warm IRAS colours and evidence for high excitation through 
the [OIII] 5007A to H/3 ratio. Again, this is true for most of 
our sample. Figure 6(b) shows the extinction corrected [OIII] 
5OO7A to Hfl ratio as a function of F60Mm/F25^im. Clearly, 
most of the HBLRs tend to have higher apparent excitation 
than the sample as a whole, and the two subsets are indeed 
statistically different at the 99.9% confidence level. This is 
not due to a dependence of the excitation on AGN luminos- 
ity as there is no correlation between the [OIII] 5007A to 
H/3 ratio and Lasoot, nor is there a tendency for the more 
luminous AGN to have larger Hell 4686A to [OIII] 5007A 



line ratios as might be expected if there was such a corre- 
lation. A large [OIII] 5007A to H/3 ratio is not a necessary 
requirement for an HBLR however: NGC5995 has a ratio 
^ 6, closer to that of the non-HBLRs with standard Seyfert 
classifications in Table 2 (cf NGC5135 and NGC7130). Nor 
is a large [OIII] 5007A to H/3 ratio a clear sign of an HBLR, 
or of warm IRAS colours, since there arc galaxies in our 
sample with cool colours and large [OIII] 5007A to H/3 ratio 
and without an HBLR. 

3.3 Radio and X-Ray Properties of the Sample 

The two other wavelengths with greatest power for exam- 
ining the properties of the active core arc cm radio and 
hard (2-lOkeV) x-ray. Shorter wavelength radio emission, 
and lower energy x-ray emission, can be due to thermal emis- 
sion from star formation if present. We therefore searched 
the literature for both integrated and core radio fluxes at 
2.3GHz, as well as the x-ray fluxes in the 2-10 keV band, and 
the neutral hydrogen column densities inferred from these. 
The data we used in our study are summarised in Table 3. 
The presence of many upper limits in the following requires 
the use of the techniques of survival analysis in comparing 
samples. We have made use of the software package ASURV, 
version 1.1 (La Valley, Isobe & Feigelson 1992), which im- 
plements methods for univariate and bivariate problems, as 
outlined in Feigelson & Nelson (1985) and Isobe, Feigelson 
& Nelson (1986). 

One caveat should be made with regard to the core 
radio fluxes. Most of these values were obtained using the 
Parkos Tidbinbilla Interferometer (PTI), which linked the 
64m Parkes antenna with the 70m Tidbinbilla antenna over 
a 275 km baseline. Most were obtained in snapshot mode, so 
if there is asymmetric structure on the scales to which the 
PTI is sensitive (~ 0.1 arcseconds) the flux will be underes- 
timated if the PTI beam is not aligned with that structure. 
A comparison of the galaxies observed by both Sadler et al. 
(1995) and Thean et al. (2000) at 8.4GHz VLA show this ef- 
fect clearly, even for some systems which are only marginally 
resolved by the PTI. However, most of our sources are likely 
to have unresolved cores, rather than more complex radio 
jet structures, and the PTI observations cover most of our 
sample, unlike the available VLA observations. Thean et 
al. (2000) did report detections of some compact cores at 
8.4GHz for sources that have only upper limits in the PTI 
observations. We have not included these since it is clear 
from work such as Sadler et al. (1995) that the spectral 
index of the core varies widely, whereas the general emis- 
sion from the galaxy is reasonably well matched to an in- 
dex near —0.75. Furthermore, the 8.4GHz emission can be 
dominated by star formation. Wc note here though that at 
8.4GHz the core in 0518-2524 contributes more than 50% of 
the total flux, whereas the cores in NGC4388 and NGC5194 
contribute less than 10%. For comparison the 2.3GHz limits 
on all 3 arc less than 10%. 

There is no significant difference in the inferred core or 
total radio power between the galaxies showing evidence for 
an HBLR and those that do not. However, there is a weak, 
though not statistically significant, trend for the HBLRs to 
have slightly higher core radio power. Unfortunately the lim- 
ited amount of data available on the core radio fluxes and 
the large number of upper limits make it difficult to say 
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whether this is real. There is also weak evidence (signifi- 
cance only 15% however) for the galaxies containing HBLRs 
to have higher ratios of core radio flux to far infrared flux 
(Figure 7a), but not higher ratios of integrated radio flux to 
far infrared flux (Figure 7b) . These plots arc essentially mea- 
sures of how important star formation is to the radio and 
far infrared emission (see Helou, Soifer & Rowan-Robinson 
1985). Galaxies in which the ratio of the total radio to FIR 
flux is high are likely to be AGN dominated. Figure 7 tends 
to indicate that most, but not all, HBLRs lie in systems in 
which the AGN dominates. 

We also examined how well the radio luminosities cor- 
related with the optical and infrared data. Lasoot correlates 
very well with both the core and total radio power at better 
than the 99% confidence level. By contrast, Lfir correlates 
with the total radio power at the 99% confidence level, but 
with the core radio power at only the 95% confidence level. 
The FIR-radio correlation for all radio quiet galaxies is well 
known (eg de Jong et al. 1985). The fact that Seyferts show 
a larger scatter in this relationship when compared to star- 
bursts is also well known (Norris, Allen & Roche 1988, Sopp 
& Alexander 1991). A discussion of the relative contribu- 
tion of the AGN core to the total flux can be found in Roy 
et al. (1994, 1998) and Heisler et al. (1998), who note that 
the scatter seen in the FIR-radio correlation is largely due 
to those galaxies with powerful radio cores. Our results axe 
consistent with these findings. In particular, the tight corre- 
lation between Lasoot and the core radio power strengthens 
the case for Lasoot being a good measure of the AGN lu- 
minosity, and the poorer correlation between Lfir and the 
core radio power indicates that L_f/_r is not solely a measure 
of the AGN luminosity. 

The 2-lOkeV x-ray band is important since it allows a 
clean measure of the intrinsic AGN flux, and many nearby 
AGN have been observed spectroscopically in that band by 
ASCA. For most of the galaxies in our sample the active core 
will be completely obscured below 2keV. The x-ray data can 
be used to derive the neutral hydrogen column densities to 
the central source. For sources with Nh < 10^'*cm~^, the 
values are derived from the observed photoelectric cut-off. 
Sources with Nh > 10^'' cm^^, often described as Compton 
thick, are obscured at lOkeV, but may show reflection dom- 
inated spectra. The 'observed' extinction in these sources is 
therefore low since none of the direct flux actually reaches 
us. The inferred high column density is determined from 
other properties, such as the ratio of the observed 2-lOkeV 
flux with extinction corrected [OIII] 5007A flux, which in- 
creases as the column density decreases, and the 6.4keV Fe 
Ka line equivalent width, which increases as the scattered 
component becomes more dominant. In a few cases, Comp- 
ton thick sources have been observed with the BeppoSAX 
satellite, which can detect emission in the 10-lOOkeV band as 
well. This allows more stringent limits to be placed in these 
cases (eg. Matt et al. 2000). The actual values as taken from 
the literature are given in Table 3, along with the original 
references. 

For two of the objects in Table 3 we have used our own 
reductions of archived ASCA data. We only used the GIS2 
and GIS3 data, and combined and grouped these data so 
that after rebinning there were at least 20 counts in each 
bin. We fitted absorbed power-laws plus a Gaussian (to rep- 
resent the FeKa line) to the data in the 2-lOkeV range, 



analogous to the procedure in Bassani et al. (1999). This is 
actually a good fit to the observed spectrum of NGC5995, 
which is relatively bright and well detected by ASCA. The 
equivalent width of the best fit line is 240ljgj5cV, and the 
photon index, F = I.S^q ]^^. The rather low obscuring col- 
umn seen is therefore intrinsic, and not due to us only seeing 
the source through scattered radiation. The ASCA data for 
IC3639 are of relatively low signal-to-noise, and allowed only 
a crude estimate of the equivalent width, of 4200l42o"cV, if 
the photon index was fixed at 2, and the line centre fixed 
at 6.34keV. However, this is consistent with the object be- 
ing Compton thick, as also found by Risaliti, Maiolino and 
Salvati (1999). We use their value of Nh, derived from Bep- 
poSAX data, since they determine a larger lower limit. 

In Figure 8(a) we plot Wasoot against Nh- There is a 
clear separation between the HBLRs and the non-HBLRs, 
in the sense that only the galaxies with the largest value 
of Wa5007 at any given obscuration show HBLRs. Figure 
8(b) shows how the sample falls in the LAsooT-A^i^ plane in- 
stead. Although it is less clear cut, it is generally true that 
the galaxies with the largest values of Lasoot at any given 
Nh show HBLRs (1925-7245 is the obvious exception here). 
Risaliti et al. (2000) reported am upper limit for the flux from 
0019—7926. We have compared this limit with other data by 
ratioing with the extinction corrected [QUI] fiux. This pro- 
vides a crude measure of the obscuring column (eg Bassani 
et al. 1999). From this we conclude that Nh > lO^^cm"^. 
We did not plot this result on Figure 8 given the uncertain- 
ties involved, but we note that this limit is consistent with 
the results from the other non-HBLRs (though see also the 
discussion in the Appendix). 

At some level the relatively smooth delineation between 
HBLRs and non-HBLRS in Figure 8(a) compared to the 
rougher trends seen in Figure 8(b) implies it is not sim- 
ply AGN luminosity, but how that luminosity dominates 
over the host galaxy that is important in the detection of 
HBLRs. It should bo noted that the separation seen in Fig- 
ure 8(a) holds true even though there are many lower lim- 
its. Essentially these data points can only move to the right, 
and not vertically. The separation between the HBLRs and 
non-HBLRs at high column density is approximately a hori- 
zontal line in this plot, so the HBLRs and non-HBLRs from 
our sample can never mix in this diagram. Of course many 
of the non-HBLRs have not been observed in the hard x- 
ray band. However, all bar three of these galaxies have ob- 
served Wa5007 <40A, and therefore are consistent with the 
trends already shown in Figure 8(a) regardless of what their 
x-ray properties might be. The exceptions axe Mrk 1361, 
NGC 5929 and NGC 7592 which have Wasoot ~60A. Hard 
x-ray observations of these galaxies would be valuable in de- 
termining whether the trends seen hold for the entire sample, 
as would spectropolarimetry of NGC 7592. 



Finally, we note there is no correlation between Nh and 
IRAS colour, in agreement with Alexander (2001). This con- 
firms the overall suspicion that IRAS colour is not directly 
linked to obscuration. 
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4 DISCUSSION 
4.1 Analysis 

It is worth stressing what our spcctropolarimctric data ac- 
tuaUy reveal before considering how wo can explain these 
observations. Clearly, the HBLRs, by definition, show ev- 
idence of an obscured (hidden) BLR, consistent with the 
predictions of the unified model. By contrast, we cannot say 
that the non-HBLRs do not contain a BLR. All we can say is 
that for the most part we do not detect scattered broad Ha 
from such a region at the level we would predict (see Section 
3.1). For most of the non-HBLRs the signal-to-noise in our 
data is sufficient that we can strongly rule our such emission 
at the predicted level. The conclusion must be that any scat- 
tered emission from the BLR, if it exists in these galaxies, 
appears at a considerably lower level than our prediction on 
the basis of the behaviour of the known HBLRs. For some 
of the non-HBLRs wo also know the BLR docs exist from 
hard x-ray observations. Although not relevant to the dis- 
cussion of why we see an optical HBLR in some galaxies, 
this does provide additional confirmation that the unified 
model is correct at some level, and indicates that the non 
detection of a HBLR is not because there is no BLR to sec. 
Clearly it would be useful to have x-ray spectra of the entire 
sample. This should now be possible, even for the very faint 
Compton thick sources, using Chandra and XMM-Netwon. 

It may also be true that with sufficient signal-to-noise 
we would discover broad Ha in spectropolarimetry in all 
Seyfert 2s. Clearly, our data shows that for many of the 
non-HBLRs the scattered flux must be significantly below 
our predicted values. Therefore only data which can actually 
probe an order of magnitude deeper in the scattered flux is 
probably of use. It may be that instrumental limitations on 
the accuracy that can be achieved in the error in the polar- 
isation will actually prevent such observations from being 
achieved even with larger telescopes. The discussion below 
allows the possibility that every galaxy contains an HBLR, 
since our basic assumption in analysing these data is that 
all Seyfert 2s do contain a BLR. If we cam actually detect 
such emission in the future, the debate would then centre on 
whether there are correlations between the detected scat- 
tered broad Ha flux and factors such as AGN luminosity 
rather than simple detectability. For the moment, however, 
we are limited to seeking an explanation as to why some 
galaxies show obvious HBLRs and others show nothing. 

There is a relatively simple explanation for all our data, 
even though our results may appear rather complex at first 
sight. First we note the clear message from Figure 3 is that 
HBLRs do appear like reddened Seyfert Is. This is evidence 
that the AGN largely dominates Ijfir in these systems at 
least. The radio data tend to support this conclusion. By 
contrast the non-HBLRs that lie along the starburst/LINER 
track clearly cannot have dominant AGN cores. Therefore 
the clear split in colour between the bulk of the non-HBLRs 
and the HBLRs is due to the relative luminosities of the 
AGN and the host galaxy as first proposed by Alexander 
(2001). The division between the two groups in the [GUI] 
5007A to H/3 ratio seen in Figure 6(b) is largely due to 
the same cause. The galaxies with prominent starbursts will 
have enhanced H/3 relative to [QUI] leading to lower ratios 
on average in the non-HBLR population. Another plausible 
reason for lower [GUI] 5007A to H/3 ratios in some of the 



higher luminosity non-HBLRs is stratification in the NLR. 
If the higher excitation lines arise nearer the nucleus, as is 
often observed for Seyferts, more obscured objects will natu- 
rally have lower [GUI] 5007A to H/3 ratios, and lower Wxsoo? 
as well. This mechanism may explain the behaviour of some 
of the Compton thick non-HBLRs for example. However, 
both the IRAS colours and the optical line ratios are essen- 
tially secondary indicators of the underlying reason as to 
why we do or do not see HBLRs. 

The dominant reason we see an HBLR in any galaxy 
is revealed when we study the galaxies at all wavelengths 
at which the AGN dominates. The HBLRs are (i) more lu- 
mmous [GUI] 5007A sources, and in particular have more 
prominent Seyfert characteristics in the optical with large 
equivalent width line emission; (ii) have more of their mid 
infrared flux arising from the region around the AGN; (iii) 
on average have brighter radio cores; (iv) have less obscura- 
tion than non-HBLRs of similar luminosity. The first throe of 
these points indicate that there are two factors present; first 
how intrinsically luminous the AGN is, and second whether 
this dominates the luminosity compared to the rest of the 
emission from the host galaxy. Figure 8(a) indicates the sec- 
ond of these factors is as important as the first. Point (iv) 
indicates that orientation is still important (at least if we as- 
sume the bulk of the obscuration arises in a torus). Clearly, 
for a galaxy that is tilted further from our line of sight, there 
will be greater extinction towards the NLR, and hence the 
observed scattered light component will be fainter. This is 
certainly consistent with the results for those non-HBLRs 
which have good signal-to-noise data, where, as noted above, 
the limit on any broad Ha present is well below the expected 
value from previously known HBLRs. 

The most obvious explanation as to why the AGN lu- 
minosity and its contribution to the bolometric luminosity 
are so vital is that more luminous AGN cores support larger 
extended scattering regions. Clearly, an extended region is 
required in those HBLRs where the obscuring column to the 
BLR is > lO^^cm^^. This region of high obscuration must 
itself be compact (cf the discussion on NGC1068 in the Ap- 
pendix and in Risaliti et al. 1999). Our data cannot constrain 
the size of the scattering region however, though imaging 
polarimetry can. Infrared imaging polarimetry of NGC1068 
and the considerably less luminous Circinus Galaxy (Pack- 
ham et al. 1997, Lumsden et al. 1999, Ruiz et al. 2000) 
clearly shows that the latter has a smaller scattering region. 
Other observable indicators also tend to favour such a corre- 
lation between the size of the scattering region and the AGN 
luminosity (see Section 4.2). If this is correct then we predict 
that the size of the scattering region in the less luminous, 
non-HBLR, Compton thin galaxies must be small. 

The alternative possibility is that we are simply limited 
by signal-to-noise in our polarimetry as claimed by Alexan- 
der (2001), and perhaps indicated by our analysis in Section 
3.1 for some of the non-HBLRs. We made an additional test 
to determine whether signal-to-noise is a major factor in 
HBLR detectability by examining the x-ray data, since the 
broad Ha flux should scale with the x-ray flux. We compared 
the x-ray fluxes in the Compton thin galaxies to test whether 
they showed a difference between the HBLR and non-HBLR 
systems. There is a clear distinction between the two groups 
in terms of luminosity, with the HBLRs being more lumi- 
nous as they are at other wavebands, but none in terms of 
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detected flux. This suggests again that luminosity is the key 
factor rather than simple signal-to-noise. It may be the case 
that those galaxies with a high current star formation rate 
also have more molecular gas near the nucleus, making it 
easier to hide the NLR. Further infrared imaging polarimc- 
try, which has the benefit of penetrating the obscuration 
whilst minimising any contribution to the polarization from 
the host galaxy, may allow us to isolate the scattering sites 
even in the less luminous or more heavily obscured systems 
and determine if it is the size of the scattering region or 
simply low signal-to-noise that best explains our results. 

Better spatial and spectral resolution would also be a 
benefit in future spcctropolarimctric surveys. Improved spa- 
tial resolution will allow the scattering region to be isolated, 
reducing the diluting contribution from the host galaxy, as 
appears to be the case for Circinus. This could be partic- 
ularly important for those cases where there is significant 
star formation. Young stars arc a likely candidate for many 
galaxies showing evidence for the so called second featureless 
contirmurn (Tran 1995b, and see also the discussion in Gon- 
zalez Delgado et al. 1998). These stars may also help mask 
any polarization signal, since they dilute the direct flux com- 
ponent from the AGN, resulting in a lower net polarization. 
Improved spectral resolution will make it easier to remove 
the contribution of polarized narrow lines, as shown graph- 
ically by Young et al. (1996a) for the case of 0518-2524. 

The presence of star formation is not in itself an ar- 
gument against orientation dependent obscuration in the 
cooler galaxies however. At least for AGN of equivalent lu- 
minosity, galaxies with significant star formation in which 
the core is obscured are more likely to be characterised as 
having starburst activity than those where there is little ob- 
scuration. This may actually be the case for at least some 
of our sample, particularly for galaxies with approximately 
equal inferred AGN luminosity such as NGC5135, NGC7130, 
NGC5995 and 0518 2524, aU of which show purely Seyfert 
indicators in the optical spectra. We know there must be 
a scattering region of some kind in the Compton thick ob- 
jects NGC5135 and NGC7130, since the hard x-ray contin- 
uum from these objects is reflected light from the Seyfert 1 
core. Again, wo should be able to test the extent, and the 
extinction to this scattering region using infrared imaging 
polarimetry. 



widths of [QUI] 5OO7A are rarely reported, it is clear from 
the published spectra in Moran et al. that these galaxies all 
have WxBooT ^ lOOA, consistent with the trends we find. The 
same is also true of the high polarization HBLRs considered 
by Awaki et al. (2000), and those HBLRs found by Young et 
al. (1996a). In fact the only exception to the trend shown in 
Figure 8(a) that we could find in the literature was for Circi- 
nus, which has Waboot ~ 30A, but Nh = 4.3 x 10^® cm"^ 
(Matt et al. 1999). It may be that the proximity of Circinus 
(so that wc have higher quality data, with better spatial res- 
olution) is the major factor in enabling us to see the HBLR. 

The case for an extended scattering region in the high- 
est luminosity sources is also made, in part, by Awaki et 
al. (2000). They show that almost all of the highly polar- 
ized HBLRs found by Miller and Goodrich (1990) and Tran, 
Miller and Kay (1992) are Compton thick, and hence that 
we see the x-ray emission as well as the optical through scat- 
tered light. This implies the scattering region must project 
outside the area of highest obscuration in all these sources. 
Their analysis indicates the x-ray emission can arise up to 
30pc away from the actual core, considerably larger than 
the putative Ipc scale height of a torus. Indeed it is possible 
that these most luminous Seyfert 2s only appear as such be- 
cause they are largely edge-on. There is certainly a dearth 
of objects in the top left hand area of Figure 8(a), and of 
all the objects considered by Awaki et al. only Was 49b lies 
in this region. Perhaps such objects are in fact classified as 
reddened Seyfert Is, or 1.5s, instead of Seyfert 2s. Lastly, we 
note that the split in x-ray behaviour reported by Awaki et 
al. between the highly polarized HBLRs from the Miller and 
Goodrich and Tran, Miller and Kay studies and their sam- 
ple of 'normal' Seyferts is in fact primarily a split between 
those systems which are Compton thick and show only rc- 
fiected emission from the core, and those in which the core 
can be seen through the obscuring material. The objects 
they report as non-HBLRs actually contain several known 
HBLRs (including four in the sample reported here), and 
therefore their comparison of the cause of the HBLR/non- 
HBLR split should be viewed instead as the split between 
Compton thick and Compton thin sources. Therefore the 
case they make that HBLRs must all be edge on is only 
strictly applicable to the highly polarized sources. 



4.2 Comparison with other work 

A detailed discussion of other observations from the litera- 
ture of some of the galaxies in our sample is presented in 
the Appendix. 

Kay and Moran (1998) and Moran et al. (2000) re- 
cently reported several new HBLRs. Examination of the 
IRAS fiuxes of these galaxies show that they too have 
warm colours (with the exception of NGC3081 for which 
no IRAS data exist). The same is true of all other previ- 
ously known HBLRs. Only three of the galaxies studied by 
Moran et al. have been observed in the hard x-ray band, in- 
cluding NGC3081, which has Nh = 6600t;^22 x lO^Ocm-^. 
Of the other two, NGC2273 has Nh > lO^^cm"^ and the 
coolest colour with F6()pm/F25Mm ~ 4.4, and NGC4507 has 
Nh = 2920±230xl0^"cm-^ and Feo^m/Fgs^m = 3.1 (where 
we have taken the neutral hydrogen column densities from 
Risaliti, Maiolino and Salvati 1999). Although the equivalent 



5 CONCLUSIONS 

We have completed a survey of the spectropolarimetric prop- 
erties of a far infrared selected sample of Seyfert 2s. Our 
sample spans a wide range of AGN core luminosity, and ap- 
pears representative of the local Seyfert 2 population. We 
conclude the following from our data: 

• The key factors determining the visibility of an 
HBLR are (i) AGN luminosity, and how much the AGN 
dominates the bolomctric luminosity of the host galaxy and 
(ii) the level of obscuration to the AGN. There is strong evi- 
dence that these are the only substantially important factors 
(cf Figure 8a). 

• We find a similar fraction of HBLRs in our sample 
as in the optically selected sample of Moran et al. (2000) of 
29%. It is likely that higher quality observations, or those 
at longer wavelengths may help to increase this fraction. 
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• There appears to be a split in the properties of our 
sample between galaxies which have a dominant host galaxy 
component and those which have a dominant AGN com- 
ponent as determined by their optical spectral properties. 
There is support for the conclusions of Alexander (2001) 
that at least the former group have overall properties that 
reflect the the host galaxy rather than the AGN. There is 
also support for the opposite point of view for the AGN dom- 
inated sources, that their overall properties are largely de- 
termined by the AGN luminosity and its obscuration along 
our line of sight. 

• Our data are still consistent with the unified model, 
though perhaps not in its simplest incarnation. For the 
subset of data for which both spectropolarimctry and x- 
ray spectroscopy exist, all show either scattered broad Ha, 
transmitted Seyfert 1 x-ray spectra or scattered hard x-ray 
continua, consistent with a central Seyfert 1 core. However, 
our results do seem to contradict the impression that has 
grown recently that only obscuration is important in deter- 
mining the observational characteristics of these moderate 
luminosity AGN (see, eg, Veron-Cetty and Veron 2000). 
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APPENDIX: COMMENTS ON INDIVIDUAL 

SOURCES 

Many of the gala^cies discussed in this paper have been the 
subject of other detailed individual studies. It is worth high- 
lighting some aspects of these studies since they illuminate 
the possible underlying physical conditions. 
0019-7926: As noted in Section 2.5 both de Grijp et al. 
(1992) and Vader et al. (1993) present spectra of this ob- 
ject showing clearer Seyfert characteristics. From the data 
in Vader et al. wc estimate that Wasoot ~ 200A, and the 
intrinsic line luminosity is 6 times larger than that re- 
ported here. This would move 0019-7926 in Figure 8 from 
the non-HBLR region into the HBLR one. We therefore sus- 
pect that new spectropolarimetric observations may detect 
an HBLR in this system. If true, this would increase the 
distinction between the HBLRs and non-HBLRs in terms of 
La6007, but not significantly change our other conclusions. 
NGC1068: This is the best studied of all the HBLRs cur- 
rently known. We now know that the likely orientation of the 
parsec scale torus in this galaxy is nearly edge on from direct 
radio imaging (Gallimore, Baum and O'Dea 1997). Unfor- 
tunately, the structure of the inner regions of NGC1068 is 
actually rather complex, since there is a small parsec scale 
radio structure aligned with the axis of this torus that is also 
visible in [OIII] images and imaging polarimctry from HST 
(Muxlow et al. 1996, Capetti et al. 1995, Capetti, Axon and 
Macchetto 1997), but the larger scale radio structure and 
visible ionisation cones arc rotated through ~ 30° relative 
to these (Capetti, Axon and Macchetto 1997). 
This complexity may however help to explain the variant re- 
quirements of the x-ray and the optical/infrared data. The 
BeppoSAX data reported by Guainazzi et al. (2000) show 
evidence for considerable variability in the hard x-ray emis- 
sion from NGC1068. They use this to constrain the loca- 
tion of the scattering medium at ~lpc from the nucleus. 
Since the Fe Ka line is fed by resonance scattering from 
a largely neutral medium, it further implies the obscuring 
torus is also ~lpc from the nucleus. Indeed, as has previ- 
ously been noted (Risaliti et al. 1999), unless the bulk of the 
obscuration lies at small radii, the implied mass of molecu- 
lar gas would exceed the observed virial mass in NGC1068 
and other nearby AGN. However, there is also clearly scat- 
tering on larger scales as detected both in the near infrared 
by Packham et al. (1997) and Lumsden et al. (1999), and 
in the ultraviolet regime by Capetti et al. (1995). Young 
et al. (1996b) found that a large scale torus (~ 200pc) is 
required to explain the infrared data, whilst still requir- 
ing the inner scattering radius to be ~ Ipc to explain the 
optical spectropolaximetry (Young et al. 1995). The solu- 
tion to this apparent dichotomy lies in the fact that there 
arc clearly two levels of obscuring source in NGC1068. The 
mid infrared polarimetry of Lumsden et al. (1999) requires 
a compact optically thick (Av ~ 100) region around the 
BLR, with a more diffuse larger scale molecular cloud with 



Plv = 30 hiding that region and the inner scattering cones. 
This latter region is also revealed in maps of the molecu- 
lar gas near the nucleus (Schinnerer, Eckart, and Tacconi 
1999). Clearly this also implies however that the reflected 
x-ray component should itself be partially absorbed. In the 
case of NGC1068, if our estimate of the visual extinction is 
correct, and the conversion to a column density is ~ 1/10 
that of the Galactic ISM (which is typical of AGN), then 
the implied Nh ~ 5 x lO^'^cm"^. This is sufficiently small 
to remain unnoticed in the existing data. NGC1068 should 
be taken as a cautionary example that local conditions can 
play a large role in the final global spectral properties of a 
galaxy. 

Finally, we note that the circumnuclcar starburst that lies at 
a radius of ~ Ikpc from the nucleus may actually dominate 
the far infrared flux from NGC1068 (Telesco et al. 1984). 
This suggests that the warmth of the infrared colours is in 
fact due to the dominance of the AGN component at shorter 
wavelengths. 

0518-2524: The higher resolution optical spectrum we used 
in determining the line ratios in Table 2 (which will be dis- 
cussed in greater detail in a paper in preparation) also shows 
that 0518-2524 has a strong Balmcr absorption line spec- 
trum in the blue. This is probably also the reason why many 
previous papers fail to detect the H/3 emission present. There 
is therefore clear evidence for a recent starburst episode in 
0518-2524. In addition, Clavel et al. (2000) show that there 
is relatively strong PAH emission in the mid infrared spec- 
trum of this object, which is usually taken as an indication 
of the 10/Ltm emission being due in significant part to star 
formation. It is a moot point as to whether the large far 
infrared luminosity is due to any star formation or the ob- 
vious AGN activity. The relatively low Wa5oo7 and failure 
to detect a compact core at longer radio wavelengths when 
a prominent core exists at shorter wavelengths may indicate 
the former however. It is also notable that this galaxy has 
the coolest mid-far infrared colours of any of the HBLRs 
presented here, even though the inferred extinction to the 
core from the x-ray data is small, which also suggests there 
may be a significant contribution to Ijfir from something 
other than the AGN. 

NGC4388: This galaxy is almost edge-on, so any obscuration 
hiding the BLR may actually be due in part to the host 
galaxy rather than a nuclear torus. Weak broad Ha has been 
detected in direct light in off nuclear positions by Shields 
and Filippenko (1996). They find that this is best explained 
as evidence for an HBLR, even though the direct evidence 
from polarimetry is actually marginal (Young et al. 1996a). 
The combination of both observations however gives greater 
confidence in the classification of this galaxy as an HBLR. 
IC3639: This galaxy has been studied in detail in the op- 
tical and ultraviolet by Gonzalez Delgado et al. (1998) and 
Gonzalez Delgado, Heckman and Leitherer (2001). They find 
evidence for a recent burst of star formation in the nucleus 
which contributes an equal amount to the AGN component 
to the overall bolometric luminosity. The starburst does not 
dominate the optical continuum however, since an older pop- 
ulation contributes a larger fraction of the light. 
NGC5135 and NGC7130: Detailed optical and ultraviolet 
studies of these galaxies were carried out by Gonzalez Del- 
gado et al. (1998) and Gonzalez Delgado, Heckman and Lei- 
therer (2001). They found evidence that the optical contin- 
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uum was dominated by a young starburst. As with IC3639, 
the starburst has a similar bolometric luminosity to the 
AGN. 

NGC7582: This galaxy showed brief evidence of broad Ha 
in direct flux (Arctxaga ct al. 1999), probably due to a low- 
ering of the obscuring column density, since there is also ev- 
idence for considerable hard x-ray variability in this galaxy 
(Schachter et al. 1998; Xue et al. 1998). Turner et al. (2000) 
summarise the evidence for a 'patchy' torus in this source, 
and the alternative possibility that the broad permitted line 
arose in a supernova as suggested by Aretxaga ct al. Turner 
et al. find that the best fit to their 2-100 kcV BcppoSAX 
data requires an obscuring column of 1.4 x lO'^'^cm^'^ which 
completely covers the source, plus a Compton thick absorber 
with Nh 1.6 X lO^^'cm"^ which only covers 60% of it. 
Because the Compton thick component does not cover the 
entire source, 'holes' in the thinner component could lead to 
unobscured sight lines towards the BLR. The optical classifi- 
cation (Table 2) is a reflection of fact that there is significant 
ongoing star formation in this galaxy. 
The galaxies without spectropolarimetry: Four of the galax- 
ies in our initial sample do not have spectropolarimetry. 
There is little available data for two of these (NGC5427 
and NGC5899) but they would appear to be typical of the 
other galaxies which have cool IRAS colours in our sam- 
ple. NGC7592 (also known as Mrk928) has been the sub- 
ject of considerable study. It consists of an interacting pair: 
the eastern component is classed as a pure starburst (Lons- 
dale, Lonsdale & Smith 1992). The western component is 
classed as a possible AGN (see Table 2), though clearly has 
characteristics of a transition or composite object. Veron, 
Goncalves, & Veron-Cetty (1997) discuss this aspect of the 
galaxy in more detail. It seems unlikely that any of these 
three gala^xies will show evidence for an HBLR. The final 
'unobserved' galaxy is 0031-2142. This system is a known 
luminous x-ray source (sec Moran, Halpern & Helfand 1996), 
but has the appearance of a starburst galaxy from its op- 
tical spectrum. Only the presence of weak broad Ha led 
to its classification as a Seyfert. The recent x-ray study of 
Georgantopoulos (2000) has confirmed its identification as 
a marginally obscured Seyfert. We have shown where this 
galaxy falls on our Afjf-WA5007 diagnostic plot. On the basis 
of this evidence it is possible that despite its cool colours this 
galaxy may show scattered broad Ha. Although the optical 
and FIR emission is dominated by the starburst component, 
the obscuration is low. In this sense it might be thought of 
as analogous to a system such as NGC5995, but with added 
star formation present. Spectropolarimetry of this source is 
certainly desirable. 
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Table 1. Complete IRAS selected sample for spectropolarimetry observations. Redshift data are taken from Strauss et al. (1992). IRAS 
fluxes arc derived from the IRAS BGS survey (Soifcr ct al. 1989; Sanders ct al. 1995) except for 0019-7926, NGC5995 and NGC5899 
which are from Strauss et al. (1990), and Mrk334, 0031-2142, 0425-0440, Mrkl361 and 1408+1347 which were taken from FSC. Galaxies 
which have a * in the flag column do not have spectropolarimetry. 
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Table 2. Optical properties of the sample. The second column indicates whether we detected a HBLR in our spectropolarimetric data. 
A n? indicates that no HBLR was detected, but that we cannot absolutely rule out the possibility that one exists with the current 
signal-to-noise of our data, as discussed in Section 3.1. A ? indicates the signal-to-noise is likely to be insufficient to determine whether 
an HBLR is present or not. A blank in this column indicates the object was not observed. The third column gives the reference for the 
optical spectropolarimetry. The codes are as follows: a - this work; b - Ruiz et al. (1994); c - Antonucci & Miller (1985); d - Inglis 
et al. (1995); e - Young et al. (1996); f - Inglis et al. (1993); g - Miller & Goodrich (1990). The other columns give the properties as 
derived from the direct fluxed spectra. The fourth column is the extinction, derived assuming the intrinsic Ha/H/3 ratio is 3.1. Columns 
5—9 give the extinction corrected ratios (as log 10 (ratio)) for (5) -'00^"^ jO^^^ ^ (g) esoaAJOl],^ 6j84A_Jnii]^ ^.g-^ 6717-H6731A [SII] 

^^^^Jlsll Column 10 gives the equivalent width of the 5007A [GUI] line in —A, and column 11 the extinction corrected flux in the 
5007A [GUI] line in units of 10" erg s ^ cm ^. Columns 12—15 give the derived classification following Veilleux & Osterbrock (1987) 
from (12) t*^™ ! versus ; (13) -^^^ versus ■'g^; (14) -^^^ versus (15) presence of Hell emission. The codes are as follows: 

I — intermediate classification; L — LINER; N — no detection; S — Seyfert. Finally, the reference for the source of the spectral data is 

given (which in some cases is not the same as the source of the spectropolarimetric data). Where more than one letter is given, the first 
refers to the source of the line ratios, and the second to the source of the equivalent width of the [GUI] 5007A line. The references are 
as follows: (1) Osterbrock & Martel (1993) - Seyfert classification is based on the probable detection of [FeX] by Osterbrock & Martel, 
which would not be expected in a typical LINER; (2) Dahari (1985); (3) this work or other unpublished data of our own; (4) Coziol et 
al. (1993); (5) Moran, Halpern & Helfand (1996) classify this as a Seyfert 1.8 on the basis of weak broad Ha; (6) Kewley et al. 2000; (7) 
Ho, Filippenko & Sargent (1997); (8) Young et al. (1996); (9) De Robertis, Hutchings & Pitt (1988); (10) Keel et al. (1985); (11) Stauffer 
(1982); (12) Golina, Lipari & Macchetto (1991); (13) Colina, Sparks & Macchetto (1991); (14) Lonsdale, Lonsdale & Smith (1992); (15) 
Kewley et al. (2001); (16) Veilleux et al. (1995). 
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Table 3. Data taken from the literature on the compactness of the lO/xm emission, C, compact radio flux, integrated radio flux (both 
in mjy, and measured at 2.3GHz), observed hard x-ray flux (in units of lO^^^ergs cm"^ and measured in the 2-lOkeV band) and 
inferred neutral hydrogen column density (in units of lO^^em^'^) from the x-ray data. Where an entry is left blank no data exists in the 
literature for that particular item for that galaxy. The 10/im compactness parameter, C, is either taken from Giuricin, Mardirossian & 
Mczzcti (1995), or derived from data presented in Maiolino et al. (1995), with the exception of NGC1068 where the flux used is taken 
from Lumsden et al. (1999). Radio data are taken from the following references: (a) Roy et al. (1998); (b) Roy et al. (1994); (e) Rush, 
Malkan & Edelson (1996); (d) Heisler et al. (1998); (e) Sadler et al. (1995); (f) Bicay et al. (1995); (g) Condon et al. (1990); (h) Morganti 
et al. (1999); (i) Su et al. (1996); (j) Condon et al. (1998); (k) Roy & Norris (1997); (1) Bransford et al. (1998); (m) Condon, Anderson & 
Broderick (1995); (n) Sopp & Alexander (1992); (o) Ulvestad & Wilson (1984). Where more than one reference is given, and a eorc flux 
is present, the core flux comes from the first listed, and the total flux from the second. If more than one reference is given, and there is 
no core flux, we have had to extrapolate the total flux from more than one source. Where the reference is marked with a *, the total flux 
has been extrapolated to 2.3GHz assuming a spectral index of -0.75. The total flux is interpolated from published data (usually at 1.5 
and 5GHz) where the reference is marked with a t. The same is true for the core flux in the case of NGC5929. X-ray data are taken from: 
(1) Risaliti et al. (2000); (2) Georgantopoulos (2000); (3) Matt et al. (1997); (4) Kii et al. (1996); (5) Isasawa et al. (1997); (6) Risaliti 
et al. (1999); (7) this work; (8) Ueno (1997); (9) Turner et al. (1997b); (10) Terashima et al. (1998); (11) Pappa, Georgantopoulos and 
Stewart (2000); (12) Guainazzi et al. (1998); (13) Malaguti et al. (1998). We have made use of the catalogue presented by Bassani et al. 
(1999) in compiling these x-ray data. 
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Figure 1: IRAS colour-colour plots of galaxies from the survey of Kcwlcy ot al. (2000). Seyfert 1 and 2s are shown in the 
left-hand panels, and starburst and LINERs (which have similar colours since LINERs tend to be dominated by the Galactic 
component as far as their IRAS colours are concerned) in the right-hand panels. The solid line is the reddening line for Seyfert 
galaxies derived by Dopita et al. (1998). The dashed line is the locus of starburst galaxies found by the same authors, and 
the dot-dashed line represents the maximal extent of the region occupied by galaxies of mixed excitation. 
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Wavelength Angstroms 

Figure 2: Spcctropolarimotry of NGC5995. The broad Ha is clearly evident in the polarized flux. Weak broad Ha is also 
present in the direct spectrum, leading to a classification for this galaxy as a Seyfert 1.9. 
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Figure 3: IRAS colour-colour plots as in Figure 1, but for the sample of Seyferts considered here. The crosses represent those 

galaxies in which an HBLR has been detected, the solid circles are those galaxies without a detection and the triangles are 
those galaxies for which we do not have spectropolarimetry. Note how the HBLRs tend to congregate on the Seyfert reddening 
line, in the region of warmer IRAS colours. Some of the Seyferts considered clearly do not have colours dominated by the 
AGN from this plot, since they have colours similar to the starburst /LINER population in Figure 1. 
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Figure 4: (a) Lasoot and (b) Lfia as a function of IRAS mid-far infrared colour. There is a weak correlation present in (a) 
but not in (b), suggesting that the Lasoot, but not Lf/a, is a factor in the observed colour. The symbols are as in Figure 3. 
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Figure 5: The correlation between the [OIII] 5007A line equivalent width and the extinction corrected [OIII] 5007A luminosity. 
The symbols are as in Figure 3. 
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Figure 6: The observed correlation between (a) the equivalent width of [OIII] 5007A, and (b) the extinction corrected [OIII] 
5007A to H/3 ratio, and the IRAS mid-far infrared colour. On average the HBLRs have larger absolute equivalent widths, 
and higher values of the optical line ratio. The symbols in (a) are as in Figure 3. 





Figure 7: The ratio of (a) core radio flux (essentially an AGN property), and (b) the integrated radio flux, and the 60/im 
flux. These are both measures of how luminous the central AGN is relative to any star formation present. The first shows a 
trend for HBLRs to have higher values of the ratio. Since the actual radio luminosities between the HBLRs and non-HBLRs 
are not significantly different, this shows that the 60/im flux must be larger in the non-HBLRs. The trend is not as obvious in 
(b) especially if the radio galaxy IC5063 is removed from consideration, showing the integrated radio luminosity is determined 
in part at least by the total activity present in the galaxy and not just the core. The symbols are as in Figure 3. 
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Figure 8: The data in our sample plotted as a function of (a) AGN to host galaxy luminosity and obscuration and (b) 
AGN luminosity and obscuration. All of the non-HBLRs lie at lower AGN luminosities than the HBLRs at equivalent column 
densities. This shows how luminosity and obscuration are both important in determining our ability to see an HBLR. The 
symbols are as in Figure 3. 
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